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The 5*-wave effective range parameters of the neutron-deuteron (rid) scattering are de- 
rived in the Faddeev formalism, using a nonlocal Gaussian potential based on the quark- 
model baryon-baryon interaction fss2. The spin-doublet low-energy eigenphase shift is suf- 
ficiently attractive to reproduce predictions by the AV18 plus Urbana three-nucleon force, 
yielding the observed value of the scattering length 2 a n a and the correct differential cross sec- 
tions below the deuteron breakup threshold. This conclusion is consistent with the previous 
result for the triton binding energy, which is nearly reproduced by fss2 without reinforcing 
it with the three-nucleon force. 
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§1. Introduction 

Few-nucleon systems are best suited to study the underlying nucleon-nucleon 
(NN) interaction and its extension to few-nucleon forces, since many sophisticated 
techniques to solve the systems yield equivalent results that can be compared with 
ample experimental data. 1 ''—' In fact, the three-nucleon (3iV) system is already 
solved with many realistic meson-exchange potentials, yielding insufficient binding 
energies of the triton missing 0.5 to 1 MeV without the 3N forced An attempt to 
reproduce the NN and 3N data consistently is pursued by using the chiral effective 
field theory, but a complete reproduction of all the 3N data in this approach is 
still beyond awayF*™ We have applied the quark-model (QM) baryon-baryon (BB) 
inter actiorP'^ to the triton and hypertriton in the Faddeev formalism and obtained 
many interesting results jHUnj The most recent model fss2 gives a nearly correct 
binding energy of the triton with the correct root-mean-square radius, preserving 
the sufficient strength of the tensor force for the deuteron and the correct 1 Sq NN 
scattering length.® This result indicates that fss2 is sufficiently attractive in the 
2 Si/2 channel of the 3N system without the 3iV force. In this channel, the deuteron 
distortion effect related to the strong short-range repulsion of the NN interaction 
is very important. In the QM BB interaction, this part is mainly described by the 
quark-exchange nonlocal kernel of the one-gluon exchange Fermi-Breit interaction, 
which has quite different off-shell properties from the phenomenological repulsive core 
described by local potentials in the standard meson-exchange models. The nonlocal 
effect resulting from the exact antisymmetrization of six quarks is very important to 
reproduce the nearly correct triton binding energy. 

The QM BB interaction is constructed for two three-quark clusters in the frame- 
work of the resonating-group method (RGM). It is characterized not only by the 
dominant nonlocality from the interaction kernel, but also by the energy dependence 
originating from the normalization kernel. In the original evaluation of the triton and 
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hypertriton binding energies, this energy dependence is determined self-consistently 
by calculating the expectation value of the two-cluster Hamiltonian with the square 
integrable three-cluster wave function.®"^ This prescription is however not appli- 
cable to the scattering problem, since the scattering wave function is not square 
integrable. In the final form of the triton and hypertriton Faddeev calculations,^ 
the energy dependence of the RGM kernel is eliminated by a standard off-shell trans- 
formation, using the square root of the normalization kernel. —^ An extra nonlocality 
emerges from this off-shell transformation as a result of eliminating the energy de- 
pendence of the RGM kernel. It is shown in Ref. 11) that this renormalized RGM 
prescription gives a slightly less attractive effect to the triton and hypertriton bind- 
ing energies, in comparison with the previous self-consistent treatment. The 350 keV 
deficiency of the triton binding energy predicted by fss2 after the charge dependence 
correction of the NN interaction is still much smaller than 0.5 to 1 MeV given by 
the standard meson-exchange potentials. It is therefore interesting to examine the 
nonlocal effect of the QM NN interaction to the 3N scattering observables, such 
as the scattering lengths, the differential cross sections and the spin polarization, in 
this energy-independent QM NN interaction. 

The numerical calculation of the three-body scattering using the QM BB inter- 
action is very time consuming due to the complex structure of the interaction. We 
therefore construct a nonlocal Gaussian potential in the isospin basis by applying 
the Gauss-Legendre integration formula to special functions appearing in the ex- 
change RGM kernel.^ The nonlocality and the energy dependence of the QM BB 
interaction is strictly preserved in the nonlocal Gaussian potential. We find that 
the 15-point Gauss-Legendre integration formula is good enough to carry out the 
few-body calculations accurately. The NN phase shifts predicted by this potential 
are essentially the same as those by fss2 with accuracy of less than 0.1°. We will 
show that the difference in the triton binding energy between fss2 and this nonlocal 
Gaussian potential is only 15 keV. 

Here, we study the low-energy neutron-deuteron (nd) elastic scattering below 
the deuteron breakup threshold based on the formulation developed in RefJ14p. For 
this purpose, the S'-wave effective range theory in the channel-spin formalism is very 
useful. The channel spin S c is composed of the NN total angular momentum / 
and the spin 1/2 of the third nucleonl® Since the deuteron channel with 1 = 1 
only survives in the asymptotic region, the scattering amplitudes for the nd elastic 
scattering are specified by the channel spin S c = 1/2 (the spin-doublet channel) 
and S c = 3/2 (the spin-quartet channel). These two channels have quite different 
characteristics with respect to the deuteron distortion effect. Namely, in the spin- 
quartet channel the incident neutron can not penetrate deep inside of the deuteron 
due to the effect of the Pauli principle, resulting in the weak distortion effect of the 
deuteron. On the other hand, the neutron can freely approach to the deuteron in the 
spin-doublet channel, causing strong distortion effects reflecting strong sensitivity to 
details of the NN interaction. In this sense, the spin-doublet scattering length 2 a n( i 
becomes an important observable to measure if the NN interaction is appropriate 
or not. It is known for a long time that a larger triton binding energy corresponds 
to a smaller 2 a n d- This linear correlation is known as the Phillips line™ and is 
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confirmed by many theoretical calculations.^"^ Since the binding energy of the 
triton is not reproduced in the NN meson-exchange potentials, the experimental 
value, 2 a n d = 0.65 ± 0.04 fm|23 i s not reproduced either. The calculated 2 a n d is 
more than 0.9 fm if only an NN force is used. It is therefore a common practice to 
add the 3N force to reproduce the triton binding energy as well as the correct 2 a n d- A 
thorough investigation of the triton binding energy and the scattering lengths, using a 
number of meson-exchange potentials and various 3iV forces, is given in Ref.[23]). We 
can expect that the nonlocal effect of fss2 leads to a good reproduction of the doublet 
scattering length since fss2 gives a large triton binding energy close to the experiment. 
In Ref. 24), an application of the nonlocal interaction to the nd scattering lengths was 
made, using the NN interaction based on the chiral constituent quark modelP^'^ 
Since their Faddeev calculation does not treat the energy dependence of the QM AW 
interaction properly, they have obtained an insufficient triton binding energy^ 1 and 
a large doublet scattering length 2 a nc i, almost comparable to the meson-exchange 
predictions. 

The most accurate method to determine the scattering lengths is to calculate 
the zero-energy scattering amplitude directly, as carried out in Refs . fT7 |) ]i~8|) .123]) J24"|) . 
etc. In this approach, however, the zero-energy nd scattering is only examined. 
More extensive study of the low-energy nd elastic scattering can be achieved in 
terms of the effective range theoryp^ i n which a pole structure existing in the 
effective-range function for the doublet-S channel should be properly taken into 
account \y e can discuss the energy dependence of the low-energy S-wave 

phase shifts in this approach. On the other hand, this method has a problem of 
numerical inaccuracy at extremely low energies below 100 keV in the center-of-mass 
(cm) system, since the solution of the basic equation becomes very singular in the 
momentum representation. 

In this paper, we start with the nonlocal Gaussian QM NN interaction and 
eliminate the energy dependence numerically in the above mentioned renormalized 
RGM formalism. We then apply this interaction to the nd scattering and solve the 
Alt-Grassberger-Sahndhas (AGS) equation 33 ' to obtain the scattering amplitudes. 
The elastic scattering amplitudes are conveniently parameterized by the standard 
eigenphase shifts and mixing parameters defined in Ref. ll5p . The spin-doublet and 
quartet S'-wave effective range parameters: i.e., 2 a n d, { 2 r e ) n d, 4 a n d and ( 4 r e ) n d, to- 
gether with the pole parameter qQ in the doublet case, are calculated by employing 
the S'-wave single-channel effective range formula. We find reasonable agreement 
of the low-energy differential cross sections with the nd experimental data. Since 
the nd data have rather large error bars, we will also evaluate the differential cross 
sections of the pd elastic scattering by employing a simple prescription, called the 
Coulomb externally corrected approximation,® to incorporate the Coulomb ef- 
fect with some modifications to the nuclear phase shifts. We will find that the 2 S 1 / 2 
eigenphase shift predicted by our model is sufficiently attractive to reproduce the 
doublet scattering length 2 a n d and the low-energy nd and pd differential cross sec- 
tions. It is possible that the correct treatment of the Coulomb effect can reproduce 
the pd differential cross sections and the polarization observables below the deuteron 
breakup threshold without introducing the 3N force. 
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The organization of this paper is as follows. In §2.1, a brief description of 
the Faddeev formalism is given for the bound-state and nd scattering problems. 
The spin-isospin factors and rearrangement factors for the permutation operator are 
explicitly given in Appendix A. In § 2.2, we recapitulate the procedure to obtain the 
eigenphase shifts and their J-averaged central phase shifts from the solutions of the 
AGS equation. The single-channel effective-range expansion is explained in § 2.3. 
The eigenphase shifts of our model are compared in § 3.1 with those of Argonne V18 
(AV18) plus Urbana 3N force, obtained by the X-harmonics technique. The nd and 
pd differential cross sections below the deuteron breakup threshold are also discussed. 
The effective range parameters are given in § 3.2, together with the analysis of the 
S-wave contributions to the nd total cross sections. The last section is devoted to a 
summary. 

§2. Formulation 

2.1. Faddeev approach to the triton and the nd scattering 
We start with the three-body Schrodinger equation 

[E-H - V* GU - Vf GM - V 7 RGM ] V = , (2-1) 

where V^ GM denotes the energy-independent renormalized RGM kernel for which 
the detailed derivation is given in Ref.[T4"|). The subscripts a, (3 and 7 in Eq. (|2- 1|) 
specify the types of Jacobi coordinates related to the residual pair in the usual way, 
with (a, P, 7) being the cyclic permutation of (123). For the systems of three identical 
particles, the Faddeev equation for the bound state reads 

il> = G tPiP , (2-2) 



where P = P(i2)-P(i3) + -P(i3)-P(i2) is a sum of the permutation operators for the nu- 
cleon rearrangement.^ The NN £-matrix in the three-body model space is derived 
from the standard Lippmann-Schwinger equation t = v + vG^t with v = y RGM ; 
where Go = 1/{E — Hq) is the three-body Green function for the free motion. It 
is composed of the negative total energy E in the cm system and the three-body 
kinetic-energy operator Hq = h,Q + /jq. The operators ho and ho correspond to the 
kinetic energy for p% = (l/2)(fei — £2) and 93 = (l/3)(2fe3 — (fei + £2)) respectively 
when the Jacobi coordinates with 7 = 3 is chosen. The vector fcj is the individual 
momentum of particle i. The Faddeev component tp is defined through tp = GqVF, 
with the total wave function & = ip a m Eq. (|2Tp . In the partial wave expansion, 
we use the channel-spin formalism specified by (I^)S C . The relative angular momen- 
tum t between the spectator nucleon and the iViV subsystem is coupled with channel 
spin S c and makes the total angular momentum (£S C )J. The AW channel is specified 
by (As)I; t, where A, s and t are the orbital angular momentum, the spin and the 
isospin of the iViV system, respectively. The iViV isospin t is uniquely specified by A 
and s from the Pauli principle (— ) A + S +* = — 1. We further set the parity restriction 
7r = (— which is conserved for each J. The angular-spin-isospin wave functions 
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are thus denned by 

|p,g;123) = £ \p,q,l)(j\p,q; 123) , 

7 

(p,g;123| 7 ) = [^(5)[[>A(p)Xrt(l ) 2)]jxii(3)]s ] JJ ,.i ra , (2-3) 

with 7 = [^[(As)/i]5 c ] JJ Z ; (ti)|T 2 in the channel-spin representation. The partial 
wave expansion of the Faddeev equation in Eq. ()2-2p is given by 

]\/f i roo />oo pco 

h 2 K f + p 2 + (3/4) ^ Jo Jo Jo 

dp"p" 2 (p, q, j\t\p', q, j')(p', q', j'\P\p", q" , j")tp 7 "(p", q"), (2-4) 

where M = (M n + M p )/2 is the averaged nucleon mass in the isospin formalism. 
The binding energy of the triton is deduced from Eb = (-E) = (Ti 2 /M)k 2 . The 
coupled integral equation for p and q in Eq. fj2-4[) is solved in the Lanczos-Arnoldi 
method after the necessary process of discretization. See Ref.[8j) for details. The 
AW i-matrix is factorized as 

^\^^^ = ^ m ^^{^ E -W^) ■ (2 ' 5) 
The matrix element of the permutation operator is evaluated as 

(p,q,~f\P\p',q',j') = ^[ dx 6< " P x+ 2^ 9-y,Y (g, q', x) ^, y ,f^ , (2-6) 
l j_l p ^ p ^ 

with pi = p(q', q/2; x), P2 = p(q,q' /2;x) and p(a,b;x) = \/ a? + b 2 + 2abx. The 
basic rearrangement coefficients g 7j y (q, q' , x) contain the spin-isospin factors and 
the explicit expression depends on a specific type of the channel-coupling scheme, as 
given in Appendix A. 

For the nd scattering, the three-body scattering amplitudes are obtained by 
solving the AGS equation^ 

U\<f>) = G 1 P\<p) + PtGoU\4>) , (2-7) 

where \cp) = \qo,if^d) is the plane- wave channel wave function with |^) being the 
deuteron wave function. The total cm energy E in Go -1 = E + iO — Ho is expressed 
as E = E cia + Ed, where E cra = (3h 2 /4M)(/q is the cm incident energy of the neutron 
and \ed\ = —£d is the deuteron binding energy. In the AGS equation, we have in 
general two types of singularities, but the notorious moving singularity does not 
appear for the energies below the deuteron breakup threshold. The other singularity 
related to the deuteron pole of the AW t-matrix is directly incorporated to the AGS 
equation in the Noyes-Kowalski method. For the detailed procedure to overcome 
difficulties of these singularities, Ref.[M|) should be referred to. Here, we recapitulate 
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the method to derive the elastic scattering amplitudes. The singularity of the NN 
i-matrix is separated as 

t = t-icGn 1 \6)(6\Ga 1 with c = 2ir^- . (2-8) 



If we use the principal- value i-matrix t in Eq. (j2-7f) . we obtain 

= G^P^fl - ic<0|!7|0>] +PtG U\<t>) . (2-9) 
We eliminate the first term of Eq. (|2-9p by defining W as 

W = G P - P\<f>)Z- x {4\P with Z=(</>\Gq 1 P\<I>) . (2-10) 
Our basic equation is 

Q\tf>) = P|0) + WtQ\<f>) , (2-11) 
where Q|0) and -P|</>) is defined through 

G o tf|0> = Q|<M<zTO>, = mz- 1 . (2-12) 



The elastic scattering amplitudes are obtained by multiplying Eq. (|2-9p with Z 
from the left-hand side. One can easily show 

{<l>\U\<t>) = [Z- 1 - {<t>\X\<t>) + ic}- 1 with {4\X\if>) = (<f>\PtQ\<l>) ■ (2-13) 

The partial wave components of the scattering amplitude, U^, S r\ ^ s \ = {<!>£' S' c \U\<t>(.s c ) 
are defined by 

{4>qp S' c S' cz \U\4> qi \ S C S CZ ) = ^2 U(i'S' c ),(eS c ) 

i'UJ Z 

x J2( i ' m ' S 'c S 'cz\JJz)Ye mf (q f ) Y,^™ScS cz \JJz)Yt m {qi) , (2-14) 



with \qf \ = \qi\ = q^. Once the partial wave components of {4>\X\4>) is calculated, 

'(l'S' e U£S c ) 



U/p,o,s (f a\ is obtained by solving an equation, 



The coupled channel S-matrix is given by S^ s n /#gi\ = &S C ,S' C ~ ^ c ^us ) (£'S'V 
2.2. Eigenphase shifts 

In the channel-spin representation, the asymptotic channel wave function 
\4>; (£S C )JJ Z ) in the partial wave expansion is specified by (£S C )J = ( J ± 3/2, 3/2) J, 
(J T 1/2, 1/2) J, and (J =F 1/2, 3/2) J for the parity vr = (-) J=Fl / 2 . Therefore, the S'- 
matrix S^ iS . ^/g/j is a two-dimensional matrix for J = 1/2, and a three-dimensional 
matrix for J ^ 1/2. The S'-matrix can be diagonalized as 

S = U^e 2iA U , (2-16) 
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where A is the diagonal matrix of the eigenphase shifts 5f Sc . (We follow the notation 
in Ref. [T5|) . but we should note that t and S c are not good quantum numbers.) The 
matrix U can be parameterized in terms of the mixing parameters s, £ and 77,^ 

/l \ / cos£ sin£ N 

U = cose sine 10 

\0 —sine cose J \— sin^ cos^y 

cos £ cos 77 cos £ sin 77 sin £ 

— cos e sin 77 — sin e sin £ cos 77 cos e cos 77 — sin e sin £ sin 77 sin e cos £ 
sin e sin 77 — cos e sin £ cos 77 — sin e cos 7/ — cos e sin £ sin 77 cos e cos £, 

(2-17) 

We assume e = £ = for the J w = 1 /2 + and £ = 77 = for the J" = 1 /2~ . By using 
a function, 

1 _ 1 + iz , 

Arctan z = —Log , (2-18) 

2i 1 — iz 

we can calculate the mixing parameters e, £ and 77 from the following equations: 

. , U12 , a 1 / / C/'3icos77 + J732sin77 
77 = Arctan , 4 = —Arctan 



Oil V ^33 

/C/siShlT? - [/ 32 COS77\ 

e = -Arctan — : , 2-19) 

\ U21 sin 77 - U22 cos 77 J 

where Uij is the component of the matrix U. 

Before applying Eq. (|2-19[) . we should reorder the eigenstates of Eq. (|2-16p such 
that the leading term of the eigenvectors follows the decreasing order. This prescrip- 
tion corresponds to the rule to assign the eigenphases to such quantum numbers 
as having the dominant components of the eigenvectors F^l For energies above the 
breakup threshold, all the angles become complex. The largest component of the 
eigenvector is given by a real positive number. This gives an ambiguity of the sign of 
the mixing angles, which in turn comes from the phase convention of the eigenstates. 
We should note that Eq. (|2- 19[) is most convenient if the diagonal Uu components 
are close to one. This is because the mixing angles are usually very small with the 
magnitude of several degrees (except for E cm > a few MeV in the 2 Pj- 4 Pj couplings). 

Since the difference in the eigenphase shifts with the same (£S C ) but different 
J is very small, the J-averaged central phase shift 5f s is convenient to discuss the 
scattering cross sections approximately. The J-averaged central phase shift is defined 
by taking an average of all possible J states 2S c +1 £j with the weight factor (2J + 1), 

^ = (2 < + l)(k + D ^ (2J+1) ^- <2 ' 20) 

The differential cross section for the nd elastic scattering is calculated by summing 
over the final spin states and by taking an average over the initial spin states. This 



*' This rule is not applied to the case of 2 f3/2 _4 f3/2 coupling for _E cm > 40 MeV, where the 
magnitude of the mixing parameter e exceeds 45°. 
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is given by 



da 

~dn 



1 (do_ 
3 \d~Q 



) 



Sc=1/2 2fda\ Sc=3/2 

+ 3 \d~a) 



(2-21) 



where the differential cross section for each channel spin S c is calculated from the 
scattering amplitudes fi = (l/qo)e lSe sin^ through 



to compare our results with the experimental data, we have to take into account the 
Coulomb force for the pd scattering. The comparison with the pd data is desirable 
since they are abundant and more accurate than those of the nd scattering. For 
example, in Refs.l38| H 43| ) many phase shift analyses have been carried out with high 
accuracy. The exact treatment of the Coulomb force in the three-body problem is 
still a challenging taskPS , '33'S31 | ''^ Because of its long-range nature, the Coulomb 
potential is not amenable to the standard scattering theory. We therefore incorpo- 
rate the Coulomb effect to evaluate the differential cross sections of the pd elastic 
scattering, by employing a following simple prescription with some modifications to 
the nuclear phase shifts. The single channel differential cross section is given by, 



where f c (0) is the standard Coulomb scattering amplitude. The scattering am- 
plitude from the strong interaction, fP in Eq. (|2-23p . is estimated from the nd 
eigenphase shifts of fss2 by adding the difference of the corresponding quantities 
for the pd and nd scattering, which will be taken from published results for other 
NN interactions. In Eq. (|2-23|) . the partial- wave Coulomb phase shift era is given by 
an = argT(£ + 1 + in) with the Sommerfeld parameter n for the relative motion of 
the proton and the deuteron. 

2.3. Effective-range expansion for the nd scattering 

If we use the effective range theory, we can study energy dependence of the phase 
shifts reflected in the effective range r e and discuss contributions to the total cross 
sections from the S-wave components between the neutron and the deuteron. We 
first calculate the eigenphase shifts for J 71 " = l/2 + and 3/2 + states. In the case of 
J 77 = l/2 + , the 2 Si/2 and 4 -Di/2 channels are coupled. For the neutron energies below 
the deuteron breakup threshold, the phase shifts and mixing parameters are small 
except for the dominant S-wave eigenphase shift, so that the effective-range expan- 
sion formula for a single channel problem can be safely applied to this component to 
obtain the effective range parameters. This is also the case for the J w = 3/2 + state, 
where the 4 5 , 3 / 2 , 2 -C*3/2 an d 4 D 3 / 2 channels are coupled. In the quartet 5-channel, 





Here, we abbreviate the J-averaged central phase shift 5f s to b~i for each S c . In order 




(2-23) 
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Fig. 1. The cut structure used in the N/D formalism for the nd elastic scattering amplitude in the 
complex z planeP^ Here, z is defined by z = E cul /\e c i\- 

we can expand the effective-range function K(qo) = qo cot 5 in the power series of q$: 

K(qv) = -- + \r e ql + 0{qt) , (2-24) 



where a is the scattering length, r e is the effective range, and 5 = 5^ 3 in K(qo) = 

°2 

go cot 8 is the S- wave eigenphase shift for the quartet channel with go being the 
relative wave number between the neutron and the deuteron. 

In the doublet-,? channel, the effective-range function K(qo), which is the real 
part of the inverse scattering amplitude, has a pole just below the elastic thresh- 
oldP3*-E2Ii Wg parameterize the effective-range function in the doublet channel as 



, fW = J±M±^M), (2 . 25) 



21 

1 + (qo/q Q Y 



i 

where a pole parameter qQ specifies the pole position, and S = 6*i is the 5-wave 

U 2 

eigenphase shift for the doublet channel. The origin of this pole structure is studied 
by means of N/D equations EB-EZJi j n jy/D formalism, the partial wave compo- 
nents of the scattering amplitude are given by N(z)/D(z), where N(z) and D(z) are 
the analytic functions of the complex and dimensionless energy variable z defined by 
z = E cm /\eci\- The two and three-body unitarities and the cut structure in the nega- 
tive energy, which are shown in Fig. 1, yield the relationship, N(z) = —z~2lm D(z). 
The N/D equations are usually constructed by applying the Kramers-Kronig rela- 
tions to N(z) and D(z). From the solution of the N/D equation, it was found that 
this singularity are brought about by both the dominant single-nucleon exchange and 
the other effects such as the two-nucleon exchange etc. slowly varying for small z. In 
the doublet channel, the single-nucleon exchange, which is by far the longest-range 
force, is attractive. Thus, nothing prevents the other effects from influencing the 
low-energy scattering. On the other hand, the single-nucleon exchange is strongly 
repulsive in the quartet channel. Therefore, the nucleons in this channel can not 
penetrate to the region where the other attractive forces could act. This is why the 
pole structure is found only in the doublet channel. 
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§3. Results and Discussion 

3.1. Eigenphase shifts 

The direct comparison of the nd eigenphase shifts, predicted by fss2, with the re- 
sults of the modern phase shift analysis for the pd scattering is not possible because of 
the Coulomb effect. We therefore list in Table U our results for the energies of E n =l, 
2 and 3 MeV together with other theoretical predictions by the Pisa group,® which 
are calculated using Argonne V18 NN potential (AV18) 46 > and AV18+Urbana(UR) 
3A^ potentials. 47 -* Here, E n = (3/2)E cm is the neutron incident energy measured in 
the laboratory system. We have included the AW interaction up to the total angular 
momentum / max = 4 and the momentum mesh points n = n\-n2-n^ =12-6-5 in the 
notation defined in §3.1 of Ref. ll4p . The UR 3A^ potential gives a sizable effect of 
about three to four degrees only on the J n = l/2 + channel. We immediately find 
an outstanding feature in the J n = l/2 + channel. Namely, our results by fss2 are 
very similar not to the AV18 results but to the AV18+UR(3A r ) results shown in the 
parentheses. It is not surprising since fss2 gives the nearly correct binding energy 
without introducing the 3A^ force. On the other hand, the phase shift parameters of 
the J T = 3/2 + state are very similar between fss2 and AV18. In this state, the effect 
of the UR 3A" force is very small owing to the Pauli principle. The difference between 
fss2 and AV18 is less than 0.2° - 0.3°, which is comparable to the effect of the 3N 
force. For the P states, some of the eigenphase shifts show somewhat larger differ- 
ence from the AV18 and AV18+UR(3A r ) results especially at E n = 3 MeV, but still 
less than 1° difference. After all, we have found good correspondence between our 
fss2 results and the predictions by the AV18+UR(3iV) potentials. The resemblance 
seen in Table U becomes more transparent if we calculate the J-averaged central 
phase shifts defined by Eq. (|2-20p and compare them. The J-averaged central phase 
shifts can be used to evaluate the nd differential cross sections through Eqs. (|2-2ip 
and (|2-22|) . We have illustrated these in Figs. [2]- 4 with dot-dashed curves, but 
they almost overlap with the solid curves for the exact calculations. We find that 
the D-wave components give an appreciable contribution to the differential cross 
sections even in such a low energy as E n =l MeV. This is of course because of the 
D-wave component of the deuteron wave function. This analysis encourages us to 
study the pd differential cross sections by a simple approximation for the Coulomb 
effect, discussed in § 2.2. We first assume that the nuclear scattering amplitude fj? in 
Eq. (|2 23j) is equal to the nd scattering amplitude from the J-averaged central phase 
shift, fe = (e 2 *^ s c — l)/2iqo. This prescription yields fairly large overestimation of the 
differential cross sections as plotted in Figs. [2]- 4 with dotted curves. We find that a 
large effect of Coulomb modification on jf is necessary for the present low-energy pd 
scattering, which is a well-known fact claimed by many authorsP^'ESffil 1 '!^ Here, 
we use an extended version of the Coulomb externally corrected approximation's in 
which the nuclear eigenphase shifts for the pd scattering, 5f s (pd) , are calculated from 
our 5f s (nd) by adding the difference of those evaluated with another AW interaction. 
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Table I. The nd eigenphase shifts and mixing parameters (in degrees), obtained from the model 
fss2. The maximum angular momentum for the NN system, J ma x = 4, and the momentum mesh 
points n = 12-6-5 are used. The corresponding parameters calculated by the Pisa group from 
the AV18 potential models are also listed for comparison.^ The parameters in the parentheses 
are predictions by the AV18+UR(37V) potentials. 



Model 


fss2 


AV18 


fss2 


AV18 


fss2 


AV18 


Sn(MeV) 


1.0 


1.0 


2.0 


2.0 


3.0 


3.0 




-0.978 


-0.980 


-2.52 


-2.53 


-3.82 


-3.85 


— 


(-0.976) 


— 


(-2.52) 


— 


(-3.84) 


2 5'i/2 


-14.8 


-18.1 


-24.2 


-28.3 


-30.8 


-35.3 


— 


(-14.3) 


— 


(-24.0) 


— 


(-30.8) 


»7l/2+ 


1.29 


0.928 


1.48 


1.08 


1.55 


1.12 


— 


(1.39) 


— 


(1.47) 


— 


(1.45) 


2 Pl/2 


—4.12 


—4.13 


—6.55 


—6.57 


—7.43 


—7.49 




i A 1 Q\ 

(-4.13) 




( — 0.08 ) 




(-7.50) 


4 Pl/2 


11.8 


12.0 


19.3 


19.9 


23.6 


24.2 




(12.1) 




(20.1) 




(24.5) 


£l/2- 


3.50 


3.47 


5.02 


4.98 


6.76 


6.68 


— 


(3.53) 


— 


(5.07) 


— 


(6.82) 


4 'S'3/2 


—46.6 


—46.7 


—60.5 


—60.8 


—69.6 


—69.9 




\— 40.0 J 




{ fin i\ 
(-60.7) 




y— oy. / ) 


2 P>3/2 


0.564 


0.564 


1.50 


1.51 


2.34 


2.36 




(U.004 J 




(1.51) 




{Z.00 ) 


i D , j,/2 


— 1.05 


— 1.05 


—2.71 


—2.72 


—4.12 


—4.14 




(-1.05) 




(-2.71) 




( A 1 A \ 

(-4.14) 


£ 3/2+ 


0.603 


0.621 


0.688 


0.686 


0.763 


0.747 




(U.OZo ) 




^U.ooo ) 




i C\ 7K /I \ 
(U. 1 t>4J 


£3/2+ 


0.516 


0.511 


0.957 


0.948 


1.36 


1.35 




(0.514) 




(0.948) 




(1.35) 


V3/2+ 


-0.106 


-0.107 


-0.232 


-0.231 


-0.363 


-0.363 


— 


(-0.105) 


— 


(-0.228) 


— 


(-0.356) 


4 TP 

F3/2 


0.122 


0.121 


0.491 


0.488 


0.919 


0.920 




(0.121) 




(0.489) 




(0.921) 


2 D 
P3/2 


—4.06 


—4.08 


—6.38 


—6.41 


—7.10 


—7.18 




( A fl£^ 
^ — 4.UO ) 




(-OA6) 




( 7 9D\ 
{— 1 .ZD) 


4 P 3 /2 


13.7 


13.9 


21.9 


22.3 


25.5 


26.0 




(14.0) 




(22.3) 




(26.0) 


£3/2- 


-1.28 


-1.24 


-1.93 


-1.86 


-2.72 


-2.62 




(-1.27) 




(-1.89) 




(-2.66) 


?3/2- 


-0.196 


-0.177 


-0.334 


-0.262 


-0.427 


-0.265 




(-0.177) 




(-0.259) 




(-0.256) 


'73/2- 


-1.04 


-1.00 


-2.17 


-2.17 


-3.57 


-3.52 




(-1.04) 




(-2.16) 




(-3.53) 


4 G 5 / 2 


-0.015 


-0.015 


-0.091 


-0.090 


-0.206 


-0.206 


2 D 5/2 


0.560 


0.559 


1.49 


1.49 


2.31 


2.33 


4 D 5/ 2 


-1.11 


-1.11 


-2.90 


-2.90 


-4.44 


-4.46 


£ 5/2+ 


-0.263 


-0.277 


-0.291 


-0.297 


-0.312 


-0.315 


?5/2+ 


-0.233 


-0.272 


-0.491 


-0.494 


-0.716 


-0.701 


V5/2+ 


-0.659 


-0.821 


-1.42 


-1.49 


-2.02 


-2.04 
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Table I. — continued 



Model 


fss2 


AV18 


fss2 


AV18 


fss2 


AV18 


£„(MeV) 


1.0 


1.0 


2.0 


2.0 


3.0 


3.0 




13.1 


13.2 


21.4 


21.7 


25.8 


26.0 


— 


(13.2) 


— 


(21.8) 


— 


(26.3) 


2 F 5/2 


-0.063 


-0.063 


-0.251 


-0.251 


-0.465 


-0.466 


F 5 /2 


0.127 


0.127 


0.514 


0.510 


0.947 


0.951 


£5/2 — 


0.399 


0.447 


0.479 


0.472 


0.514 


0.538 


^5/2— 


0.384 


0.390 


0.690 


0.684 


0.938 


0.926 


%/2- 


-0.117 


-0.123 


-0.239 


-0.239 


-0.343 


-0.334 


4 -C*7 /2 


-1.03 


-1.03 


-2.66 


-2.67 


-4.04 


-4.06 


2 

G7/2 


0.0076 


0.0075 


0.047 


0.047 


0.108 


0.107 


4 

G7/2 


-0.015 


-0.015 


-0.094 


-0.095 


-0.215 


-0.214 


£7/2+ 


0.230 


0.325 


0.353 


0.368 


0.355 


0.355 


?7/2+ 


0.362 


0.427 


0.781 


0.798 


1.16 


1.14 


^7/2+ 


-0.106 


-0.143 


-0.287 


-0.299 


-0.459 


-0.459 


^7/2 


-0.062 


-0.062 


-0.248 


-0.248 


-0.458 


-0.460 


4 F 7/2 


0.132 


0.132 


0.533 


0.532 


1.00 


1.00 


£7/2- 


-0.207 


0.238 


-0.244 


-0.236 


-0.256 


-0.232 


2 Gg/2 


0.075 


0.0075 


0.047 


0.046 


0.106 


0.105 


G9/2 


-0.015 


-0.016 


-0.097 


-0.097 


-0.223 


-0.223 


£9/2+ 


-0.125 


-0.183 


-0.189 


-0.199 


-0.173 


-0.176 


4 F g/2 


0.124 


0.124 


0.497 


0.496 


0.921 


0.922 




-0.015 


-0.015 


-0.092 


-0.091 


-0.208 


-0.208 



Table II. The nd J-averaged central phase shifts (in degrees), calculated from the eigenphase shifts 



in Table I through Eq. (|2-20|). 


. The others 


are the same 


as Table I. 






Model 


fss2 


AV18 


fss2 


AV18 


fss2 


AV18 


E n (MeV) 


1.0 


1.0 


2.0 


2.0 


3.0 


3.0 


■ 2 s 


-14.8 


-18.1 


-24.2 


-28.3 


-30.8 


-35.3 






(-14.3) 




(-24.0) 




(-30.8) 




-4.08 


-4.10 


-6.44 


-6.46 


-7.21 


-7.28 






(-4.10) 




(-6.48) 




(-7.30) 


2 D 


0.562 


0.562 


1.49 


1.50 


2.32 


2.34 


2 F 


-0.0626 


-0.0624 


-0.249 


-0.249 


-0.462 


-0.463 


4 S 


-46.6 


-46.7 


-60.5 


-60.8 


-69.6 


-69.9 






(-46.6) 




(-60.7) 




(-69.7) 


4 P 


13.1 


13.2 


21.2 


21.6 


25.3 


25.7 






(13.3) 




(21.7) 




(25.9) 


4 D 


-1.05 


-1.05 


-2.73 


-2.74 


-4.15 


-4.18 


4 F 


0.127 


0.127 


0.510 


0.508 


0.948 


0.950 
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do/dr2(mb/sr) E N =1 MeV 



600 



400 



200 



100 



~\ solid nd fss2 




exp. 

pd: Ko69 
nd: EI62 



60 120 

9c.m.(deg) 



180 



Fig. 2. fss2 predictions to the nd and pd differential cross sections obtained from various prescrip- 
tions for the phase shifts at En = 1 MeV: the exact nd calculation (solid curve), the J-averaged 
central phase shifts (dot-dashed curve), the Coulomb externally corrected approximation with 
5 — 8(nd) (dotted curve), and the Coulomb modified nuclear phase shifts in Table III (dashed 
curve). The experimental data are taken from Ref. 48) for E162 (nd with errorbars) and Ref. l49j> 
for Ko69 (pd with filled circles) . 



do/dn(mb/sr) En= 2 MeV da/dfi(mb/sr) E N = 3 MeV 




e c . m .(deg) 6 cm (deg) 

Fig. 3. The same as Fig. 2, but for £^=2 MeV. Fig. 4. The same as Fig. 2., but for E N =3 

MeV. The experimental data are taken from 
Ref. l50fl for Sc83 (nd with errorbars) and 
Ref. l51fl for Sa94 (pd with empty circles). 
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Table III. The Coulomb modified J-averaged central phase shifts (in degrees), S N = Sfs c (nd) ± At , 
used for the calculations of the pd differential cross sections in Figs. 2-4. Here, Sf Sc (nd) are 
given in Table Ull and Ai = [Sfs a (pd) ~ <^s c ( nc 0]Avi8 are evaluated from Ref.l42p. The results of 
the phase shift analysis (PSA) at En = 3 MeV are also shown for comparison. 



£iv(MeV) 




1.0 




2.0 




3.0 


3.0 (PSA) 


3.0 (PSA) 




At 




A t 


5 N 


A t 


5 N 


Ref.O 


Ref.HH 


■ 2 s 


4.9 


-9.9 


4.1 


-20.1 


3.1 


-27.7 


-24.85 ±0.23 


-24.87 ±0.28 


A S 


9.7 


-36.9 


7.8 


-52.7 


6.8 


-62.8 


-63.80 ±0.11 


-63.95 ± 0.28 


4p 


-3.1 


10.0 


-2.6 


18.6 


-2.0 


23.3 


23.86 ±0.01 


23.37 ±0.11 



The Coulomb correction can be evaluated since the nd and pd eigenphase shifts for 
E N = 1, 2 and 3 MeV by the AV18 potential (or the AV18±UR(3iV) potentials) 
are both given in Tables 1 and 2 of Ref.l42p. respectively. The J-averaged central 
phase shifts calculated from these values imply that the Coulomb modification is 
significant (more than 1°) only for the 2 S, 4 S and 4 P channels. The difference Ai = 
[dfs (pd) ~ fifs ( nc 0]AVi8 listed in Table ITTT1 is obtained in this way. The pd differential 
cross sections calculated with these Coulomb modified nuclear scattering amplitude 
ff are plotted in Figs. [2]- 4 with dashed curves. We find excellent agreement with 
the experimental data for the pd differential cross sections with the aid of the AV18 
Coulomb effect. The overestimation in the case of the original Coulomb externally 
corrected approximation is improved mainly by 7° - 10° modification of the A S phase 
shift to the attractive direction (see Table III). These analyses imply that the correct 
treatment of the Coulomb force could reproduce the pd experimental data for the 
differential cross sections without reinforcing fss2 with the 3N force. 

3.2. The nd effective range parameters 

Since the good correspondence between fss2 and AV18±UR(3./V) are found in the 
eigenphase shifts, we can expect that the S'-wave effective range parameters for the 
spin-doublet and quartet channels should also be reproduced by fss2 without intro- 
ducing the 3A^ force. These are determined from Eq. ()2-24p for the quartet channel 
and Eq. f|2-25j) for the doublet channel. The S'-wave effective range parameters are 
obtained by using Schlessinger's point method^ 1 (a type of Pade approximation) 
to the effective-range function K(qo) = q$ cot<5. This method is convenient to ap- 
proximate a function with a pole such as the doublet effective-range function and 
to take into account the contributions of the higher order terms in Eqs. (|2-24p and 
(|2-25p in a natural way. It is very hard in the momentum representation to maintain 
sufficient accuracy of the eigenphase shift for the extremely low energies if we use 
the realistic deuteron wave function including the D-w&ve component. We therefore 
use the sample points of energies between E cm = 200 keV and 2 MeV shown in 
Figs. 5 (the doublet channel) and 6 (the quartet channel) unless the denominator 
in the Schlessinger's point method hits zeros. In order to obtain the eigenphase 
shifts with accuracy of less than 0.01°, we need to take fine mesh points of p for the 
AW relative motion. We use a set of mesh points n=6-10-5 and the AW partial 
waves up to G-wave (/ max = 4). A typical example of fitting to K(q<f) is shown 
in Figs. 5 and 6. We reconfirm a pole structure existing in the doublet channel at 
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q cot S (fm 1 ) 




0.5 1 1.5 2 

Eam. (MeV) 



Fig. 5. The effective-range function K(qo) = 
qo cot S for the doublet S'-state as a func- 
tion of -E cm . The curve shows the rational 
approximation made by the Schlessinger's 
point methodpl 



q cot S (fm ) 
-0.08, ■ 




Fig. 6. The same as Fig. 5, but K(qo) = 
go cot 5 for the quartet 5-state. 



E Q = -(3/i 2 ^/4M) ~ -150 keV. 

Table HVl lists the model-space dependence of the triton binding energy and effec- 
tive range parameters. As to the almost converged triton binding energy £b( 3 H) = 
8.311 MeV with 

/max — 6, the small difference from the original result of fss2 (8.326 
MeV with J max = 6^P comes mainly from the nonlocal Gaussian approximation to 
the interaction kernel, adopted in this study. Since we have three NN pairs in the 
triton, this 15 keV difference in the triton binding energy is consistent with the 4 
keV difference in the deuteron binding energy, which is the difference of 2.2250 MeV 
(the original result of fss2^P and 2.2206 MeV (the result of the nonlocal Gaussian 
potential based on fss2). We find that the spin-quartet scattering length 4 a nc i is 
quite insensitive to the model space adopted. This insensitivity is related to the 
small distortion effect of the deuteron due to the Pauli repulsion of the nd interac- 
tion in this channel, which is a kinematical constraint from the spin-isospin quantum 
numbers. The system is therefore independent of the details of the NN interaction. 
On the other hand, the spin-doublet scattering length 2 a n d is subject to a strong 
channel coupling effect as seen in Table HVl We find that the simplest five-channel 
(S + D) calculation, incorporating the 3 Si + 3 Di and 1 5'o NN channels only, yields 
a value very close to the converged one. However, this is quite accidental and the 
well converged value is achieved after many partial waves, up to the G-wave of NN 
interaction at least, are included. The values of \E Q \/( 2 a nd ) ~ 220 keV/fm are al- 
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Table IV. The triton binding energy _E_g( 3 H) and S-wave effective range parameters, 2 a„d, ( 2 r e ) n d, 
4 a„d and ( 4 r e ) n( ;, predicted by fss2 for various model spaces with the maximum angular mo- 
mentum of the AW interaction (7 ma x) included. The nonlocal Gaussian potential with 15-point 
quadrature is used for fss2. In "S + D" , the 3 Si ± 3 Di and 1 So NN channels are only included. 
The pole energy for the doublet channel, Eq, is also shown. The calculated deuteron binding 
energy is 2.2206 MeV. The momentum mesh points with n=6-10-5 are used. 



^max 


£; S ( 3 H) (MeV) 


2 a nd (fm) 


( 2 r e ) nd (fm) 


E Q (keV) 


4 a nd (fm) 


( 4 r e )nd (fm) 


S + D 


8.247 


0.65 


-149 


-147 


6.30 


1.84 


1 


7.948 


0.94 


-102 


-207 


6.31 


1.85 


2 


8.213 


0.72 


-133 


-163 


6.30 


1.84 


3 


8.298 


0.67 


-146 


-151 


6.30 


1.84 


4 


8.307 


0.66 


-148 


-148 


6.30 


1.84 



Table V. Comparison of the nd scattering lengths, predicted by using fss2 (J max = 4), with other 
models. For the fss2 results, the charge dependence of the NN force is neglected. The heading 
AW implies the calculation using only the AW force, and AW±TM99 the calculation including 
the Tucson-Melbourne 99 (TM99) 27r-exchange 3iV forceP)' 1 ^ The results by CD-Bonn 2000, 
AV18 and Nijm I for the AW force are taken from Ref.l23|l (7 ma x = 5). The experimental values 
are taken from Ref, 125p . The values of 4 a n d are insensitive to the 3 A force. 



model 


£s( 3 H) (MeV) 




2 a nd (fm) 


4 a nd (fm) 




NN 


AW+TM99 


NN 


AW+TM99 


AW(+TM99) 


fss2 


8.307 




0.66 




6.30 


CD-Bonn 2000 


8.005 


8.482 


0.925 


0.569 


6.347 


AV18 


7.628 


8.482 


1.248 


0.587 


6.346 


Nijm I 


7.742 


8.485 


1.158 


0.594 


6.342 


exp 




8.482 




0.65 ±0.04 


6.35 ±0.02 



most independent of the model space. This linear correlation is already suggested in 
Ref. 132]) for the separable potentials. A strong correlation between Eb^H) and 2 a n d 
is also apparent in Table iTVl 

In TablelVl we compare £'b( 3 H) and 2 a n d by fss2 with other calculations using 
meson-exchange potentials including the 3N forceps Here, fss2 does not incorporate 
the charge dependence of the NN interaction, but the other calculations in TablelVl 
include this effect. As in £^( 3 H), it should give an appreciable influence to 2 a n d- We 
will estimate the maximum shift of 2 a n( i by simply assuming the same correlation 
as Phillips line for the triton binding energy. For fss2, the slope of the Phillips line, 
—0.686 fm/MeV, yields 0.13 fm for the charge dependence effect of the triton binding 
energy 190 keVF^J After the charge dependence correction, 2 a n d for fss2 would turn 
out to be 2 a n d ~ 0.76 - 0.80 fm. TablelVl shows that the effect of the 3N force is 
more important than the charge dependence of the NN force. When NN meson- 
exchange potentials are only used, 2 a n d is more than 0.9 fm. The experimental 
values for 2 a n d and I?b( 3 H) are reproduced only when the 3A^ force is included. The 
model fss2 almost reproduces £'b( 3 H) and 2 a n d simultaneously without introducing 
the 3N force. We should keep in mind that the mechanism to reproduce 2 a n d is 
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Table VI. The quartet (S c = 3/2) and doublet (S c = 1/2) S-wave contributions to the total cross 
sections a tot , calculated from the effective range parameters for fss2 with / max = 4 in Table IIVI 
The calculated total cross sections otot are taken from Ref. I14|) . 





MeV 


1 MeV 


2 MeV 


3 MeV 


4 S (mb) 


3325 (99.5%) 


2055 (72%) 


1467 (59%) 


1128 (54%) 


2 S (mb) 


18.2 (0.54%) 


122.4 (4.3%) 


158.4 (6.3%) 


165.5 (7.9%) 


5-wave (mb) 


3343 


2177 (77%) 


1625 (66%) 


1294 (62%) 


<Ttot (mb) 


3185 (0.15 MeV) 


2833 


2480 


2104 


exp. (mb) 


3120 ± 


2893.6 ± 18. 2 s3 


2550.6 ± 11. i 53 


2158.0 ± 7.2 s3 




(0.07 MeV) 


2854 ± 39^ 


2537 ± lcP 


2240 ± 9(P* 






3110 ± 20CP 1 


2600 ± 8(P* 


2160 ± 86^* 



quite different from that in the spin-quartet case. The large positive value for 4 aw 
is related to the Pauli repulsion for the loosely-bound deuteron cluster. On the 
other hand, the value of 2 a nc i is entirely from the dynamical origin related to the 
fairly large triton binding energy £^ xp ( 3 H) = 8.482 MeV and the existence of the 
pole structure in the effective-range function just below the elastic threshold. In 
this special situation, it is natural that more attractive nd interaction, afforded by 
the 3N force in the meson-exchange potentials or by fss2 without the 3N force, can 
reproduce a smaller value for 2 a n d, which is the content of the Phillips line. 

In Table IVT1 we show the S- wave contributions from the quartet (S c = 3/2) and 
doublet (S c = 1/2) channels to the total cross sections, calculated from the effective 
range parameters for fss2 with / max = 4 in Table IIVI We find that the total cross 
sections are dominated by the 5-wave contribution, which is more than 60% even for 
E n = 3 MeV. Furthermore, the quartet state is far more important than the doublet 
state due to the small values of \qo cot 5\, even considering the statistical weight factor 
(2S C + 1). As the energy increases, the contribution from the doublet state becomes 
appreciable owing to avoiding the pole structure just below the elastic threshold, 
but is still less than 10% at E n =3 MeV. This implies a very special situation that 
an extra attraction to the 2 Si/ 2 state by the 3iV force is unimportant to reproduce 
the differential cross sections of the low-energy nd scattering, and they are mainly 
determined by the magnitude of the repulsive 4 5 3 / 2 eigenphase shift. 

§4. Summary 

Motivated by the success of the QM baryon-baryon interaction fss2 reproducing 
the triton binding energy almost correctly without the 3iV forceps we have ex- 
tended the Faddeev calculation to the low-energy nd scattering by employing a new 
algorithm 14 ' to solve the Alt-Grassberger-Sandahs equation.^ 1 The QM AW inter- 
action, formulated in the two-cluster RGM, has rich contents of nonlocality that the 
standard meson-exchange potentials do not possessP In addition to the dominant 
nonlocality from the RGM interaction kernel, an extra nonlocality emerges from the 
off-shell transformation to eliminate the energy dependence connected to the nor- 
malization kernelp^ 1 which is sometimes neglected in similar works .EU-ESJ T/o reduce 
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the computation time for three-body calculations, we have developed and used the 
nonlocal Gaussian potential constructed from the model fss2p21 | The nonlocality and 
the energy dependence of fss2 is, however, strictly preserved in this potential model, 
resulting in an almost the same amount of the triton binding energy with only 15 
keV less. 

In this paper, we have examined the effective range parameters of the nd scat- 
tering and the low-energy differential cross sections below the deuteron breakup 
threshold. The elastic scattering amplitudes in the channel-spin representation are 
parameterized by the eigenphase shifts and mixing parameters,^ 1 from which the 
5-wave effective range parameters are derived by employing the single-channel ef- 
fective range formula. We have reconfirmed that the improved single-channel effec- 
tive range formula should be used for the channel-spin doublet (S c = 1/2) state, 
to incorporate the pole structure of the effective-range function, existing just be- 
low the elastic threshold.^ 1- ® 1 The predicted effective range parameters by fss2 

are: 2 a nd = 0.66 fm, ( 2 r e ) nd 150 fm, E Q = -(3ft 2 <?|/4M) 150 keV, and 

4 GW = 6.30 fm, ( 4 r e ) n d = 1.84 fm without the charge dependece of the NN inter- 
action. After the charge dependence correction, 2 a nd for fss2 would turn out to be 
2 «nd ~ 0.76 - 0.80 fm. It is found that the almost same 2 a n d value is accidentally 
obtained in the restricted model space involving only the 3 Si + 3 D\ and 1 5'o AW 
interaction. However, the deuteron distortion effect in the spin doublet channel is 
so strong that the sufficient partial waves, up to the G-wave at least, are necessary 
to obtain the converged result. On the other hand, the positive value A a ndi ~ 6.3 fm, 
implying the repulsive nature of the nd interaction in the spin quartet (S c = 3/2) 
channel, is quite insensitive to the expansion of the model space, owing to the kine- 
matical constraint by the effect of the Pauli principle. 

A detailed comparison of the nd eigenphase shifts predicted by fss2 and by the 
AV18 plus the UR 3A potential^ shows a prominent resemblance especially for 
the 2 S*i/2 state. The UR 3A potential gives a sizable effect only on the J 71 " = l/2 + 
channel, while very small effects on the 4 S*3/2 state and the other partial waves. 
We find reasonable agreement with the nd experimental data for the low-energy 
differential cross sections. Since the pd data are more precise than the nd data, we 
have also examined the differential cross sections of the pd elastic scattering. We have 
employed a simple prescription adding the Coulomb amplitude to the nd scattering 
amplitudes with the factor e l ( a e+o- e i) ^ w \ 1 [ c \ i [ s ca lled the Coulomb externally corrected 
approximation.® 'S3 The assumption using the nd scattering amplitude with no 
modification gives too large differential cross sections for E„ = 1 - 3 MeV, which 
implies that the Coulomb modification of the nuclear phase shifts is very important 
in this low-energy region. The Coulomb modified nuclear phase shifts are evaluated 
by adding the major difference of the nd and pd eigenphase shifts predicted by the 
AV18 potential.^ The modification is carried out with respect to the J-averaged 
central phase shifts only for the 2 S, 4 S* and 4 P channels. We have found agreement 
with the experimental data for the pd differential cross sections at E p = 1, 2 and 3 
MeV. Since the S3/2 contribution is dominant in the differential cross sections below 
the deuteron breakup threshold, the 4 £3/2 eigenphase shift of the nd (and probably 
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pd) elastic scattering is properly predicted by fss2. 

Based on these analyses, we can conclude that the spin-doublet low-energy eigen- 
phase shift predicted by fss2 is sufficiently attractive to reproduce predictions of the 
AV18 plus Urbana 3N force, yielding the observed value of the doublet scatter- 
ing length and the correct nd and pd differential cross sections below the deuteron 
breakup threshold. These results are in accordance with the bound-state calculation 
of the triton, in which fss2 predicts a nearly correct binding energy close to the ex- 
perimental value without introducing the 3N force It is important to examine 
if the correct treatment of the Coulomb force for fss2 can reproduce the pd differ- 
ential cross sections and the polarization observables below the deuteron breakup 
threshold. 
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Appendix A 

The spin-isospin factors and the rearrangement coefficients in Eq. \2- b}) 

In this appendix, we give the rearrangement coefficients g~ y (q, q , x) in Eq. 
(|2-6p for the permutation operator P = P(i23) + P( 2 i23)' together with the spin- 
isospin factors for the 3A^ system. The two contributions from -P(i23) = P(i2)P(23) 
and -P( 2 12 3) = P(i3)P(23) become equal owing to the Pauli principle for the particle 1 
and particle 2. We therefore only need to calculate the rearrangement coefficients 
for (— 2)P( 13 ). These are given by 
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and -Pfc(x) being Legendre polynomials. The spatial rearrangement factor tf\fy(\'£') 
is given by 

(2A + 1)!(2A' + 1)! 



^(\eum - ^(\>i>u\e) ~ 47r ^ + 1)!(2A 2 + 1)!(2A; + 1)!(2A' 2 + 1)' 
dqd& [Y (XlX2)x (q',q)Y t (q)]* LM P k (q-q') \r {K x i)x ,(q,?)Y e @) 
(2A + 1)!(2A' + 1)! 



LM 



(2A 1 )!(2A 2 )!(2A' 1 )!(2A' 2 )! 



XIX' t' ^(A 2 0£0|/0)(A / 2 0£ / 0|/ / 0)(A;0Ai0|/ / 0} 



x(feOA' 1 0|/0) 



A A 2 



A' a; 



Ai / k 



(A-3) 



Here, A = V2A + 1 etc. and A 2 = A - Ai, A' 2 = A' - X[ with Ai = - A, A' x = - A'. 
The explicit expression of the spin-isospin factors {Xn)~J depends on the channel 
specification scheme. For the LS coupling scheme 



(p,g;123|7) 



Y( M ) L (p,q) [x<rt(l»2)xi |(3) 



(A-4) 



with 7 = [{Xt)L(s\)S} JJ Z ; (t\)\T z , we do not need the LS sum in Eq. (|A-2|) . since 
the orbital angular momentum L and the total spin S are both conserved. We modify 
Eq. (|MD to 



( y \LSJ 
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where a common definition of the spin and isospin factors 
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(A-6) 



is used. In Eq. (|A-6p . the upper row (or the left-most column) corresponds to s=0 
(s'=0) and the second row (or the right-most column) corresponds to s=l (s'=l). 
For the jj-coupling scheme 

(p,g;123| 7 )= [[r A (^Xrf(l,2)]j[^(^xii( 3 )]il„ - > ( A ' 7 ) 
with 7 = [(As)J(^|)j] JJj; (t§)§ T z , we use the superposition of Eq. ()A-2jl with 
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The coefficients in the channel-spin formalism with Eq. (|2-3|) are similarly obtained 

as 
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LSJ 



st,s't' 



(A-9) 



For the practical calculations, it is convenient first to calculate the product of one 
6-j and one 9-j coefficients defined through 



A LSJ -X^(i\i+e+S c +jc^ is 
A (es c )(\ s )i - zJ* l > cJ \j 



5 S c I 

J j £ 



A s I 

I I 3 
LSJ 



(A40) 



Furthermore, the factor 2 will be taken off to cancel with the 1/2 factor over the x 
integral. Then we use for Eq. (|A-2p 



1 AiA^fc 
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